Universality in the flooding of regular islands by chaotic states 
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We investigate the structure of eigenstates in systems with a mixed phase space in terms of their 
projection onto individual regular tori. Depending on dynamical tunneling rates and the Heisenberg 
time, regular states disappear and chaotic states flood the regular tori. For a quantitative under- 
standing we introduce a random matrix model. The resulting statistical properties of eigenstates as 
a function of an effective coupling strength are in very good agreement with numerical results for a 
kicked system. We discuss the implications of these results for the applicability of the semiclassical 
eigenfunction hypothesis. 
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I. INTRODUCTION 

The classical dynamics in Hamiltonian systems shows 
a rich behaviour ranging from integrable to fully chaotic 
motion. In chaotic systems nearby trajectories separate 
exponentially in time and ergodicity implies that a typ- 
ical trajectory fills out the energy-surface in a uniform 
way. However, integrable and fully chaotic dynamics are 
exceptional [l| as typical Hamiltonian systems show a 
mixed phase space in which regions of regular motion, 
the so-called regular islands around stable periodic orbits, 
and chaotic dynamics, the so-called chaotic sea, coexist. 

For quantized Hamiltonian systems the fundamental 
questions concern the behaviour of the eigenvalues and 
the properties of eigenfunctions, especially in the semi- 
classical regime. From the semiclassical eigenfunction hy- 
pothesis US 0, IS El one expects that in the semiclas- 
sical limit the eigenstates concentrate on those regions 
in phase space which a typical orbit explores in the long- 
time limit. For integrable systems these are the invariant 
tori. In contrast, for ergodic systems almost all orbits 
fill the energy shell in a uniform way. For this situa- 
tion the semiclassical eigenfunction hypothesis is proven 
by the quantum ergodicity theorem which shows that 
almost all eigenstates become equidistributed on the en- 
ergy shell 0|- 

For systems with a mixed phase space, in the semi- 
classical limit (h — ► 0), the semiclassical eigenfunction 
hypothesis implies that the eigenstates can be classified 
as being either regular or chaotic according to the phase- 
space region on which they concentrate. This is sup- 
ported by several studies, see e.g. d, [l(J EL O EH- 
It is also possible, that the influence of a regular island 
quantum mechanically extends beyond the outermost in- 
variant curve due to partial barriers like cantori and that 
quantization conditions remain approximately applicable 
even outside of the island However, it was recently 
shown that the classification into regular and chaotic 
states does not hold when the phase space has an infinite 
volume [l4j. In this case eigenstates may completely ig- 
nore the classical phase space boundaries between regular 



and chaotic regions. 

In order to understand the behaviour of eigenstates 
away from the semiclassical limit, i.e. at finite values of 
the Planck constant h, one has to compare the size of 
phase-space structures with h. Let us consider for sim- 
plicity the case of two-dimensional area preserving maps 
and their quantizations. Regular states of an island con- 
centrate on tori which fulfill the EBK-type quantization 
condition 
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for the enclosed area [15| . This quantization rule explic- 
itly shows that regular eigenstates only appear if h/2 is 
smaller than the area A Iog of that island. 

Another consequence of finite h in systems with a 
mixed phase space is dynamical tunneling (l6| . i.e. tun- 
neling through dynamically generated barriers in phase 
space, in contrast to the usual tunneling under a poten- 
tial barrier. Dynamical tunneling couples the subspace 
spanned by the regular basis states, corresponding to the 
quantization condition |T]), with the complementary sub- 
space [17| composed of chaotic basis states. This raises 
the question whether the eigenstates of such a quantum 
system can still be called regular or chaotic. 

In Ref. [l8| it was shown that ([T|) is not a sufficient 
condition for the existence of a regular eigenstate on the 
TO-th quantized torus. In addition one has to fulfill 
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where T H ch = h/A c h is the Heisenberg time of the sur- 
rounding chaotic sea with mean level spacing A c h and j m 
is the decay rate of the m-th regular state, if the chaotic 
sea were infinite. When condition ([2]) is violated one ob- 
serves eigenstates which extend over the chaotic region 
and flood the m-th torus [la]. To distinguish them from 
the chaotic eigenstates that do not flood the torus, they 
are referred to as flooding eigenstates. For the limiting 
case of complete flooding of all tori, the corresponding 
eigenstates were called amphibious [l4[ ■ Recently, the 



2 



consequences of flooding for the transport properties in 
rough nano- wires were studied p^ |. 

The process of flooding was explained and demon- 
strated for a kicked system in Ref. [l8(. Condition ([2]) 
was obtained by scaling arguments, which cannot pro- 
vide a prefactor. Moreover, for an ensemble of systems, 
one would like to know the probability for the existence 
of a regular eigenstate. In particular, when varying the 
Heisenberg time, how broad is the transition regime dur- 
ing which this probability goes from 1 to 0? Another 
question is, how do the chaotic eigenstates turn into 
flooding eigenstates for a given torus? 

In this paper we give quantitative answers to these 
questions. We study the flooding of regular tori in terms 
of the weight of eigenstates inside the regular region 
and devise a random matrix model which allows for de- 
scribing the statistics of these weights in detail. Ran- 
dom matrix models have been very successful for obtain- 
ing quantitative predictions on eigenstates in both fully 
chaotic systems and s yste ms with a mixed phase space, 
see e.g. [& HJ HH, [HM H- For the present situation 
we propose a random matrix model which takes regular 
basis states and their coupling to the chaotic basis states 
into account. The only free parameters are the strength 
of the coupling and the ratio of the number of regular to 
the number of chaotic basis states. From this model the 
weight distribution for eigenstates is determined. 

For a kicked system we define the weight by the projec- 
tion of the eigenstates onto regular basis states localized 
on a given torus m. The distribution of the weights allows 
for studying the flooding of each torus separately. The 
resulting distributions are compared with the prediction 
of the random matrix model and, after an appropriate 
rescaling, very good agreement is observed. This agree- 
ment shows explicitly the universal features underlying 
the process of flooding, giving a precise criterion for the 
existence or non-existence of regular, chaotic, and flood- 
ing eigenstates in mixed systems. 

The text is organized as follows. In section [TT| we in- 
troduce the kicked system used for the numerical illus- 
trations, both classically (part A) and quantum mechan- 
ically (part B). In section [II] C we define the weight of an 
eigenstate by its projection onto regular basis states and 
investigate the distribution of the weights for the kicked 
system. In section IIIII we introduce the random matrix 
model and determine the corresponding weight distribu- 
tion as a function of the coupling strength. In section 
IIVI the relation between parameters of the kicked system 
and the random matrix model is derived. This allows for 
a direct comparison of the distributions. In section fVl we 
consider the fraction of regular eigenstates, both for an 
individual torus and for the entire island. In section I VII 
we briefly discuss the consequences of the random matrix 
model on the number of flooding eigenstates. A summary 
and discussion of the eigenfunction structure in generic 
systems with a mixed phase space is given in section IVlII 



II. THE KICKED SYSTEM 
A. Classical dynamics 

For a general one-dimensional kicked Hamiltonian 

oo 

H(p,q,t)=T(p) + V(q) £ S(t - n) , (3) 

n— — oo 

the dynamics is fully determined by the mapping of po- 
sition and momentum (q n ,Pn) at times t = n + + just 
after the kicks 

q n +i = q n + T'(p n ) , (4) 

Pn+i = p n -V'(q n+ i) . (5) 

Choosing the functions T'(p) and V'(q) appropriately, 
one can obtain a system with a large regular island and 
a homogeneous chaotic sea. For the system considered in 
[l8| , first introduced in [3] , one starts with the piecewise 
linear functions (see Fig. [TJd) 

t'{p) = i + Q - sp\ sign (p-[p+ 1/2J ) , (6) 
v'(q) = - rg -(l-r)L«+l/2j , (7) 




-0.5 0.0 p 0.5 -0.5 0.0 q 0.5 



FIG. 1: (color online) (a) Classical dynamics of the kicked 
system given by Eqs. @ and <(Sj . Invariant tori of the regular 
island are shown (continuous curves) and the transport to the 
right is indicated by the arrows. The dots correspond to one 
chaotic orbit. The magnification shows that the boundary of 
the island to the chaotic sea is rather sharp with only very 
small secondary islands, (b) Smoothed functions T'(p) and 
V'(q) (blue, dark lines) and discontinuous functions t'(p) and 
v'(q) (red, light lines) according to Eqs. (|6]|9]). 
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where [^J is the floor function, and s and r are two 
parameters determining the properties of the regular is- 
land and the chaotic sea. Usi ng a Gaussian smoothing 
with G e (z) — exp(— z 2 /2e 2 )/V27re 2 , one obtains analytic 
functions 



T'(p) = / Azt\p + z)G £ {z) , (8) 

J — OG 

V'(q) = / dzv'(q + z)G e (z) . (9) 



By construction, these functions have the periodicity 
properties 



T'(p + k) = T'{p) , 
V'(q + k) = V'(q)-k 



(10) 
(11) 



for any integer k. We consider p S [— 1/2, l/2[ and 
q e [—1/2, — l/2+M[ with periodic boundary conditions. 
The phase space is composed of a chain of transporting 
islands centered at (q,p) — (k, 1/4) with < k < M — 1 
that are mapped one unit cell to the right (see Fig. QJ,). 
The surrounding chaotic sea has an average drift to the 
left as the overall transport is zero [25, 26] . The fine scale 
structure at the boundary of the island to the chaotic sea 
has a very small area (see the magnification in Fig. [1^,). 
Resonances in this layer are irrelevant in the h regime 
studied here. For s = 2, r = 0.65 and e = 0.015 the 
regular island has a relative area A reg as 0.215. 



B. Quantization 

In kicked systems, the quantum evolution of a state 
after one period of time 



(12) 



is fully determined by the unitary operator, see e.g. [la . 

S3, IM SI, 11, 



U = exp 



2?ri .A / 27ri , , 
—V{q) exp -—T{p) 
n c fi J V "cff 



(13) 



Here the effective Planck's constant h e g is Planck's con- 
stant h divided by the size of one unit cell. The eigen- 
states of this operator are defined by 
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(14) 



where the eigenphase ipj is the quasienergy divided by 
fruj. In order to fulfill the periodicity of the classical dy- 
namics in p direction, the quantum states have to obey 
the quasi-periodicity condition 



(p + 1|V> =e- 2 ^(p|^) . 



(15) 



One can show that this leads to quantum states that are a 
linear combination of the discretized position states | (?,•), 
with qj = h c g(j + x P )- Additionally, imposing periodicity 
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FIG. 2: (color online) (a) Eigenphases of the kicked system 
vs Xq f° r ^efl = 1/10. The pattern of straight lines (inter- 
rupted by avoided crossings) with negative slope corresponds 
to regular eigenstates with m = and m = 1 whose Husimi 
functions are shown to the right. The other eigenstates are 
chaotic and live outside of the regular region, as can be seen 
from the Husimi representation, (b) Weights Wo and W\ of all 
eigenstates vs \q (left). Distribution P(W) of these weights 
in a log- linear representation (right). 



after M unit cells in q direction, quantum states have to 
fulfill the property 



{q + M\ij}) = e 2 * ix « (q\4>) 



(16) 



Because of the required periodicity the phase space is 
compact and the effective Planck's constant can only be 
a rational number 



Ks ■ 



M 
iV 



(17) 



We consider the case of incommensurate M and N, so 
that the quantum system is not effectively reduced to 
less than M cells. 
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The properties (TO}, (JTTJ) of T'(p) and V'(q) imply for 
their integrals 



T(p+k) = T(p) , 

V{q + k) = V (q)-kq~^ 



(18) 
(19) 



From this one finds that the propagator U is consistent 
with the periodicity conditions (fl~5|) and (fl"6|) if and only 
if 



(20) 



For given M and iV, this condition limits the possible 
values of the phase \p, while Xq remains arbitrary. Thus, 
in the basis given by the position states \qj), with < j < 
iV — 1, where iV is the dimension of the Hilbert space, the 
propagator U is represented by the finite N x N unitary 
matrix 



U, 



kl 



j N-l 

3=0 



' [V (?* )+T( Pi )+ K (q, -q k )} 



(21) 



where < k, I < N — 1 and pj = (j + Xq)/N. Finding 
the solution of (fT4")l . i.e. the eigenphases and eigenstates 
of the system, therefore reduces to the numerical diag- 
onalization of the matrix (f2"Tj) . The result is illustrated 
in Fig. [5]Ja) for h c g = 1/10, where the eigenphases are 
plotted as a function of Xq- The straight lines with neg- 
ative slope correspond to the regular eigenstates [H, [26| , 
whose Husimi distributions are shown to the right in 
Fig. Ufa). Lines with an average positive slope corre- 
spond to chaotic eigenstates. 

When the system consists of M unit cells one has M 
regular basis states localized on the m-th torus. Their 
EBK eigenphases are equispaced with a distance 1/M 



C. Projection onto regular basis states 

In order to investigate the amount of flooding we use 
the projection of the eigenstates onto regular basis states 
of the island region. For the considered kicked system 
regular basis states can be constructed from harmonic 
oscillator eigenstates, as the invariant tori are accurately 
approximated by ellipses 31| . The expression for the ra- 
th harmonic oscillator state, centered in a phase space 
point (q,p), is 



1 



{^) 1/4 H m (^(q-q) 



V2 m m! 

x ex P ( -it-M - ?) 2 + tLP^ - 9/ 2 ) 



(22) 



where H m is the Hermite polynomial of degree m. The 
complex constant a takes into account the squeezing and 



rotation of the state. From the linearized map at the 
stable fixed point of the island one finds a = (V351 — 
13i)/40. 

For a chain with M identical cells, a regular basis state 
is a linear combination of the harmonic oscillator states 
I^jPim)) centered in the fc-th island for < k < M — 1 
and properly normalized and periodized in the q and p 
directions [3l| . The subspace spanned by these M regular 
basis states is the same as the one spanned by the M 
harmonic oscillator states I^JT^m)- Therefore, we define 
the weight W m of a normalized state \^) by its projection 
onto this subspace corresponding to the m-th quantized 
torus 



w m = ]T |(*K 1/4 )I 2 



(23) 



fc=0 



By means of the weight W m for all eigenstates of 
Eq. (|2"Tj) we can study the process of flooding for each 
torus separately. This allows for a detailed analysis and 
a quantitative comparison with a random matrix model. 
Therefore this is a considerable improvement compared 
to our previous analysis [l8j |. where the weight was de- 
fined as the integral of the Husimi distribution of an 
eigenstate over the whole region of the island, which 
means that the information on individual tori is not ac- 
cessible. 

In Fig. [DJb) we show the weights Wq and W\ of all the 
eigenstates as a function of Xq- F° r Wo we observe that 
for almost all Xq the weights are essentially zero or one. 
Only at avoided crossings of regular and chaotic eigen- 
states their weights have intermediate values. For m = 1 
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FIG. 3: (color online) Distribution of Wo (Eq. (|23j0 vs system 
size M for effective Planck's constant h e g ~ 1/10. 
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FIG. 4: (color online) Distribution of Wi (Eq. (J23J)) vs system 
size M for effective Planck's constant h c s ps 1/10. 



the avoided crossings are much broader due to the larger 
coupling and the value W — 1 is not reached between 
several avoided crossings. This is also seen in the weight 
distributions shown to the right in Fig. [2fb), where the 
two peaks from the chaotic eigenstates (at W = 1) and 
from the regular eigenstates (at W = 0) are broader for 
to = 1 in comparison with to = 0. Note, that in the situa- 
tion of isolated avoided crossings the involved eigenstates 
are often referred to as hybrid states. 

The distribution of the weights W m allows for studying 
the process of flooding in a quantitative way. To violate 
condition ([2]) we need to increase the Heisenberg time, 
while keeping the tunneling rates j m constant. We can 
achieve this by choosing a sequence of rational approx- 
imants M/N of h cS = l/{d + g), with d S N and the 
golden mean g = (V5- l)/2 pa 0.618. This ensures that, 
while the system size M is increased, h c g is essentially 
kept at a fixed value, and therefore the tunneling rates 
7 m are independent of M . Simultaneously, the dimen- 
sionless Heisenberg time r H ch = 1/A c h increases linearly 
with M, 



N, 



ch 



1 



M 



(24) 



where we used A c h = 1 /iVch and N c h = N 
the number of chaotic states. Here TO max is the maxi- 
mal number of regular states in a single island accord- 
ing to the EBK quantization condition |T]), TO max = 
[A reg /h c s + 1/2J . As discussed in Ref. [Hj], T HjC h may be 
bounded, due to localization effects: For M larger than 
the localization length A the effective mean level spacing 
A ch ~ (AAWM)- 1 leads to r H , ch ~ XN ch /M pa Xh eS , 



where A is measured in multiples of a unit cell and N c h/M 
is the number of chaotic states per unit cell. For trans- 
porting islands, like in the model studied here, A ~ I/70 
is unusually large [3, HI, H3 > leading to a maximal value 
T H ,ch ~ Kff/ Jo- 
in Figs.[3]and[4]we show the distribution of Wq and W\ 
for d = 9 (giving approximants h e s = 1/10, 2/19, 3/29, 
5/48, . . .) for increasing system size M. For small system 
sizes we increased the statistics by varying the phase % q 
in the quantization, as it was shown in Fig. EJb). To 
present the results in a compact form each histogram 
is shown using a color scale. The horizontal strips for 
M = 1 in Fig. [3] and Fig. [4] correspond to the histograms 
previously shown in Fig. [2jb). 

In Fig. [3] one clearly observes for small M two separate 
peaks corresponding to chaotic eigenstates at W = 34 1 
and regular eigenstates with m = at W = 1. With 
increasing system size these regular eigenstates disappear 
while the weight Wo of the chaotic eigenstates starts to 
increase and they turn into flooding eigenstates. 

Comparing Fig. 0] for W\ with Fig. [3] for Wq one ob- 
serves a qualitatively similar behavior. The difference is 
that the regular eigenstates with to = 1 disappear for 
much smaller system size M sa 100 than the eigenstates 
with to = 0, as expected from Eq. (|2|) and their ratio of 
tunneling rates, 70/71 <C 1. 

For the largest values of M only flooding eigenstates 
are left which fully extend over the chaotic sea and the 
regular island. The flooding is complete and the N eigen- 
states are equally distributed in the Hilbert space. Pro- 
jecting them onto the M regular basis states leads to the 
average value W cq = M/N = h e s ~ 1/10, in agreement 
with the observed position of the peaks in Figs. [3] and [4] 
and the findings in Ref. [3| ■ 



III. RANDOM MATRIX MODEL 

The similarity in the behavior of the histograms in 
Fig. [3] and Fig.[3J suggests a universality in the process of 
flooding, which should allow for a random matrix mod- 
elling. Such models have been used successfully for the 
case of mixed systems to describe the level splitting in the 
context of chaos assisted tunneling, see e.g. [8L I22L I35L l36j| . 
In our case, we want to describe the statistics of eigenvec- 
tors for the situation of a chain of A^og regular islands. 
Here one has iV reg equispaced regular levels correspond- 
ing to the TO-th quantized torus and N c ^ COE distributed 
chaotic levels coupled by dynamical tunneling, see Fig. [5] 
For this situation we propose a random matrix model 
with the following block structure 



H 



1 1 


V 


v T 





(25) 



This matrix is chosen to be real symmetric because the 
kicked system under consideration obeys time reversal 
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FIG. 5: Schematical plot of the regular levels with spacing 
A rcg coupled with strength i>A c h to the COE distributed 
chaotic levels with mean spacing A c h. 



symmetry. As a consequence of the block structure, the 
free parameters of this model are the ratio N leg /N c h of 
the number of regular and chaotic basis states and the 
strength v of the coupling. 

The first block H lcg models the regular basis states 
associated with one specific torus, while for simplicity we 
neglect the regular basis states quantized on other tori. 
As discussed at the end of Sec. [TT] B, in the considered 
kicked system, the EBK eigenphases of the N TCg regular 
basis states are equispaced. To mimic this behavior we 
consider for i/ rog a diagonal matrix with elements (k + 
x)/iV rC g, k = 0, 1, . . . , N rcg — 1. The parameter \ can be 
chosen from a uniform distribution between zero and one. 
The energies lie in the interval [0,1] with fixed spacing 

A reg — 




0.0 W eq 0.5 T/T/ 1.0 



FIG. 6: (color online) Distribution of weights in the random 
matrix model vs coupling strength v. The ratio N leg /N c h 
approximates the value l/(8+g). The dashed lines at W — 0.5 
and W — 0.5W cq ~ 0.052 separate chaotic, flooding, and 
regular eigenstates. 



The block H c h models the AT c h chaotic basis states, 
where we assume 7V c h > N rcg . It is also a diagonal matrix 
whose elements {Ei} are the eigenphases of an N c ^ x 
iV c h matrix of the Circular Orthogonal Ensemble (COE). 
These energies {£"/} lie in the interval [0,1] with a uniform 
average density and show the typical level repulsion of 
chaotic systems. The mean level spacing of these basis 
states is A c h = l/Mh- Note, that a GOE matrix for 
this block would have been less convenient as it leads 
to a non-uniform density of levels according to Wigner's 
semicircle law. 

The off-diagonal block V accounts for the coupling 
between the regular and chaotic basis states. It is a 
N lcg x N c h rectangular matrix, where each element is 
a random Gaussian variable with zero mean and vari- 
ance (j)A c h) 2 . The positive parameter v is the coupling 
strength in units of the chaotic mean level spacing A c h- 
Thereby the results become asymptotically independent 
of the dimension N tot = N rcg + N c ^ of the matrix for 
fixed v and N Teg /Nck. 

We identify the regular region with the subspace 
spanned by the first -/V reg components. Therefore, for any 
normalized vector (* , ■ ■ • , ®N tBe -i, ^N tes , ■ *JV tot -i) 
we define the weight W inside the regular region as 

W= £ l^f. (26) 

For a particular realization of the ensemble through the 
numbers {E{\, Xi an d the block V, we compute the 
weights W of the eigenvectors. We take for the statis- 
tics only those eigenvectors whose eigenenergies are in 
the interval [0.1, 0.9] to avoid possible border effects. 
We determine the distribution of W by averaging over 
many different realizations. Increasing the matrix size 
iVtot for a fixed ratio N^g/N^ we find that the distribu- 
tion converges. Considering a ratio N rcg /N c h = 1/(8 + g) 
and a small coupling strength u rs 0.1 the distribution 
converges around -/V to t = 200. For v « 1 bigger matri- 
ces of A/tot ~ 1000 are necessary. For v w 10, we used 
Mot ~ 10000. The limiting distributions depend sensi- 
tively on the coupling strength v. 

In Fig. [6] we plot the distribution of W for different 
values of v. We have to distinguish between the un- 
coupled regular and chaotic basis states of our model 
and the resulting eigenstates in the presence of the cou- 
pling. The eigenstates fall into three classes: a) Regu- 
lar eigenstates (W > 0.5), which predominantly live in 
the regular subspace. The remaining states, which pre- 
dominantly live in the chaotic subspace, are divided into 
two classes, depending on the strength of their projection 
onto the regular subspace compared to the equilibrium 
value W cq — N rcg /N to t- This leads to b) flooding eigen- 
states (0.5Wo q < W < 0.5), and c) chaotic eigenstates 
(W < 0.5Weq). Note, that the constants 0.5 in these 
definitions are arbitrary. 

From the energy scales in the random matrix model, 
see Fig. [SI we expect three qualitatively different situa- 
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tions for the distribution of W: 

i) v -C 1, regular and chaotic eigenstates: In this 
regime the regular and chaotic blocks are practically de- 
coupled as the coupling wA c h is much smaller than the 
mean spacing of the chaotic basis states, vA c h <C A c h- 
Two sharp peaks are observable, one at W ss due to 
the chaotic eigenstates, and the other at W w 1 due to 
the regular eigenstates. The latter peak has a smaller 
weight as the density of regular basis states is smaller. 

ii) v w 1, chaotic and flooding eigenstates: Here the 
coupling ?;A c h is approximately of the same order as 
the mean chaotic spacing A c h- All regular basis states 
are strongly coupled to several chaotic basis states and 
none of the eigenstates is predominantly regular. On 
the other hand one has different types of eigenstates as 
uA c h < A rcg : Chaotic basis states, which are close in en- 
ergy to a regular basis state, strongly couple and thus 
turn into flooding eigenstates. In contrast, there are 
many chaotic basis states which are far away from any 
regular basis state and only couple weakly. These lead 
to chaotic eigenstates which show essentially no flooding 
(W < 0.5W cq ). 

hi) v ^> iVch/iVreg, flooding eigenstates: All chaotic ba- 
sis states are strongly coupled to the regular basis states, 
vA c h ^> A rog . The resulting eigenstates equally flood the 
regular subspace. The distribution of W gets a Gaussian 
shape with mean value W eq = N reg /N to t and a decreasing 
width. 

In the transition from situation i) to ii) the two peaks 
of P{W) near W = and W = 1 broaden and move to 
the center. The regular peak broadens faster, and at v w 
0.25 its maximum disappears. At v rs 1 practically no 
eigenstates are localized in the regular subspace. When 
moving from situation ii) to iii) the different types of 
chaotic and flooding eigenstates transform into a single 
type of flooding eigenstates with a similar weight W = 
Weq in the regular subspace. 

How do the resulting distributions depend on the ra- 
tio iVrog/iVch? First, the average of P(W) is given by 
W eq = N TCg /N tot = 1/(1 + N ch /N rcg ). Secondly, the 
regular peak in situation i) is independent of N^g/N^ 
apart from a trivial scaling of the normalization with 
N ieg /N c h. Numerically we checked that this is even true 
up to v w I for the distribution with W > 0.5 and 
iVreg/A^ch < 1/(8 + g). Decreasing N Tes /N c ^ enlarges the 
size of the transition regime between ii) and iii) . In par- 
ticular, the peak near W = should stay there up to 
larger values of v. 



IV. COMPARISON 

The distribution of weights for the random matrix 
model, Fig. [SJ shows a clear similarity to the results ob- 
tained for the kicked system, Figs. [3] and [4] In order to 
obtain a quantitative comparison one has to determine 
the relation between the coupling strength v of the ran- 
dom matrix model and the system size M of the kicked 



system. This can be deduced from Fermi's golden rule in 
dimensionless form 



7 - (2tt) 



2(V 2 ) 



(27) 



where the decay rate 7 of a regular state to a contin- 
uum of states with mean level spacing A is given by 
the variance of the coupling matrix elements V. In 
the random matrix model we have (V 2 ) = (fA c h) 2 , 
A = A c h = 1/iVch, and therefore (12"?| implies 



ch 



2tt 



(28) 



Applying this relation to the kicked system, we first 
note that the tunnelling rate j m for each torus can be de- 
termined numerically |3ll| (for recent theoretical results 
see [13, IH, [H, |H SI lU . The determination of the cor- 
rect value iV c h for the kicked system requires a detailed 
discussion: A regular basis state on the m-th torus, in the 
case where the tori m* , m* + 1, TO max — 1 are already 
flooded, will couple effectively to N — m*M states for 
h c s = M/N . A change of m* affects iV c h and therefore 
v. This dependence, however, can be neglected for the 
numerical comparison in our case: The ratio of the max- 
imal and minimal possible values of v is approximately 
V^l-MAl- Ac g ). For h eS m 1/10 and A rcg = 0.215 
this gives a difference of less than 7%. Therefore we sim- 
ply use the maximal value iV c h = N — M in the following. 

For these values of 7 and iV c h in Eq. (|28|) the m-th 
torus of the kicked system has a coupling strength 



_ y / 7m(l//t c ff- 1) r— 

2vr 



(29) 



This allows for rescaling the results of the kicked system 
shown in Figs. [3] and [4] from M to v using the values 
70 = 0.0015 and 71 = 0.030 [311 ] . The comparison with 
the results from the random matrix model is shown in 
Fig. [7] The agreement is very good for both tori over a 
wide range of coupling strengths v showing the univer- 
sality of the flooding process. For v > 5, however, the 
distribution reaches a constant width in Fig. [TJb) , while 
the variance decreases for the random matrix model, 
Fig. Efc). We attribute this discrepancy to the local- 
ization of eigenstates in the kicked system for M > 1000 
[T|. As a consequence, the effective number of chaotic 
basis states near an island saturates (see the discussion 
after Eq. I[24p). leading to an effective saturation of v. 

In Figs. [5] and [5] we compare individual histograms for 
the weights Wq for m = 0. To visualize the low values 
of the distributions we choose a logarithmic-linear repre- 
sentation in Figs. [5Ja) andj^a). For M — 144 one can 
distinguish the peak near W = 0, due to chaotic eigen- 
states, from the second peak caused by regular eigen- 
states. For M — 1597 these two peaks have merged and 
only a very small fraction of regular eigenstates is left. 
In both cases the distributions agree very well with the 
prediction of the random matrix model using v according 
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FIG. 7: (color online) Distributions of the weights (a) Wo and (b) Wi, taken from Figs. [3] and |4j with M rescaled to v according 
to (|29|) . (c) Result for the random matrix (RMT) model from Fig. [5] on the same scale for a better comparison. 




FIG. 8: (color online) Distribution of (a) Wo and (b) In Wo FIG. 9: (color online) Distribution of (a) Wo and (b) In Wo 
for fteff = 144/1385 (dark lines). Results of random matrix for h e a = 1597/15360 (dark lines). Results of random matrix 
model for v = 0.218 and N TCg /N ch = 1/8.618 (light lines). model for v = 0.726 and N Icg /N ch = 1/8.618 (light lines). 



to Eq. To resolve the peak near W = we show 

in Figs. |SJb) and^b) the distributions of In Wo- Again 
very good agreement with the predictions of the random 
matrix model is observed. 

Fig. [TU] shows the distribution of In W\ for m = 1 of 
all eigenstates for h e g — 13/125. We observe discrepan- 
cies at weights smaller than 10 -3 in comparison to the 
random matrix model. This difference can be explained 
as follows: Among all the eigenstates of the kicked sys- 



tem there are regular eigenstates localized on the torus 
m = which are not considered in the random matrix 
model for m = 1. These eigenstates have a negligible 
overlap with the regular basis states with m = 1 because 
they are practically decoupled and only influence the his- 
togram at very small weights. This is confirmed by com- 
puting the distribution, under exclusion of all eigenstates 
with Wo > 0.5. The resulting distribution matches re- 
markably well with the prediction of our random matrix 
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FIG. 10: (color online) Distribution of In Wi of all eigenstates 
(thin line) for h c s — 13/125. After excluding states with 
Wo > 0.5 (blue, dark line) much better agreement with the 
random matrix model (red, light line) is found. 



model. 



state, expressed in terms of tunneling rate and Heisen- 
berg time, can be transformed using Eqs. (f25)) and 
into the condition 



v < 



1 

2^ 



(30) 



The position of v = 1/(2%) is indicated in Fig. [TT1 and 
roughly corresponds to 93% of regular eigenstates still 
existing (by the W > 0.5 criterion). While in Ref. [18j 
condition |(5J) for the existence of regular eigenstates was 
obtained from a scaling argument which does not provide 
a prefactor, our random matrix model analysis shows 
that it is quite close to 1. 

For the transition regime 1/2-k < v < 1 this model 
shows a decreasing probability for the existence of a reg- 
ular eigenstate. For v > 1, which implies 



7m > (2tt) 2 



1 



(31) 



FRACTION OF REGULAR EIGENSTATES 



we find that almost no regular eigenstate exists on the 
m-th torus. Thus v = 1 defines a critical system size M m 
associated with each quantized torus 



A more global quantity than the individual distribu- 
tions P(W) is the fraction of regular eigenstates. This 
has been studied in Ref. [HI for the total number of reg- 
ular eigenstates as a function of the system size. With 
the projection onto individual regular basis states we are 
now able to study this fraction for each torus m sepa- 
rately. For the kicked system with M cells there are at 
most M regular eigenstates localized on the m-th torus. 
However, during the process of flooding, some of these 
eigenstates disappear. Thus, we define the fraction / reg 
of regular eigenstates on the m-th. torus as the number 
of eigenstates with weight W m > 0.5 divided by M. For 
small system sizes this fraction is averaged over several 
different phases \ q . To compare the resulting dependence 
on M for different values of m and h e g we determine the 
coupling strength v using Eq. (j2"9"j) . These results are 
shown in Fig. [TXJ 

For the random matrix model we compute / reg as the 
number of eigenstates with W > 0.5 divided by the num- 
ber of regular basis states N rcg , averaged over many re- 
alizations of the ensemble. As discussed at the end of 
section IIII1 the distribution P(W) for W > 0.5 is in- 
dependent of -/Vreg/Wch, apart from a trivial rescaling. 
Therefore the resulting curve f Te g{v) is independent of 
the ratio A/r eg /iV c i, in contrast to the individual distribu- 
tions. The agreement of the fractions determined for the 
kicked system with the random matrix curve in Fig. [TT] 
is very good. This shows that f reg (v) is a universal curve 
describing the disappearance of regular eigenstates. For 
v < 0.1 the fraction of regular eigenstates is larger than 
98%. For v > 1 the fraction of regular eigenstates is 
less than 1% and the corresponding regular torus is com- 
pletely flooded. 

The criterion ^ for the existence of a regular eigen- 



47T 2 fo cff 

7 m (l - h cS ) 



(32) 



With the knowledge about the flooding of individual 
tori we can now consider the total fraction of regular 
eigenstates. The regular tori with larger m have typi- 
cally a larger tunneling rate, 70 <§C 71 <C . . . <C 7m max -i- 
Therefore the flooding of the regular tori happens sequen- 
tially from the outside of the island as the system size 
increases, as found in The total fraction of regular 
eigenstates is defined as the number of eigenstates 
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FIG. 11: (color online) Fraction of regular states / reg vs cou- 
pling strength v for random matrix model (full line) and 
kicked system for various h c s and m (symbols), where the 
system size M is rescaled to v according to Eq. (|29p . Fraction 
of flooding eigenstates /n(u) for the random matrix model 
(dashed line) for N TCg /N c h = 1/(8 + g) showing a broader 
transition. 
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FIG. 12: Total fraction of regular states J-" rcg vs system size 
M according to the prediction Eq. (j33|> . (lines) in agreement 
with the data for the kicked system for h e g « 1/10 (circles) 
and heft ~ 1/30 (squares). The arrows indicate the critical 
system sizes M m according to Eq. (|32|l . 



N les /N c h — 1/(8 + 5). Note, that this fraction is de- 
fined as the number of flooding eigenstates divided by 
the number A^ot of all eigenstates. At v = 1 all regu- 
lar eigenstates have disappeared, however, the fraction 
of flooding eigenstates is just 70%. The remaining eigen- 
states are chaotic, which have no substantial weight in 
the regular subspace. For larger values of v they turn 
into flooding eigenstates. This roughly happens when 
each chaotic basis state is coupled to at least one regular 
basis state, i.e. when wA ch = A rog /2, see Fig. [5] This 
gives v — N ch /(2N TCg ) 4.8 which is in good agreement 
with the saturation observed in Fig. QT] This shows that 
the fraction of flooding eigenstates explicitly depends on 
the parameter Areg/Ach in contrast to the fraction of reg- 
ular states f Teg (v). 

Applying this result of the random matrix model to 
the kicked system where v = N ch /(2N TCg ) sa N/(2M), 
we find using Eqs. (|28|) and (|24|) . that the fraction of 
flooding eigenstates is saturated at fa = 1 for 



with weights W m > 0.5 for any m, divided by the to- 
tal number of eigenstates N. With Eq. (|32|) we get the 
prediction 



^reg(M) 



M 
AT 



m=0 




(33) 



where f re g(v) is the universal curve from the random ma- 
trix model. For small system sizes M < M m for all m 
the total fraction of regular eigenstates is ,F r eg(M) = 



Mm max /N 



ireg, 



as expected from the semiclassical 



eigenfunction hypothesis. Fig. [121 shows !F Icg (M) with a 
succession of plateaus and drops before each critical size 
M m . Considering that the ratio of successive M m only 
varies moderately, the overall behavior of T-[ eg is an ap- 
proximately linear decrease on a logarithmic scale in M , 
explaining the observations of Ref. [l8[ . The agreement 
of Eq. (j3"3"]) with the fraction of regular eigenstates for 
the kicked system for different h e s as seen in Fig. [12] is 
remarkably good. 

We conclude this section with the remark that due to 
the independence of /regt^) on the ratio N rcg /N c h one 
can obtain this universal curve by considering a simpler 
random matrix model, where only one regular basis state 
is coupled to an infinite number of chaotic basis states 
[35| . For this simpler model it might be possible to obtain 
analytical expressions for f Teg (v). 



VI. FRACTION OF FLOODING EIGENSTATES 

The random matrix model also allows for investigating 
the fraction of flooding eigenstates. While the regular 
eigenstates disappear with increasing coupling strength 
v, more eigenstates turn into flooding eigenstates with 
O.SWcq < W < 0.5. Fig. [11] shows the increasing frac- 
tion of these states for the random matrix model with 



7™ > 7— 



7T 



1 



(34) 



Note, that this prefactor increases in the semiclassial 
limit leading to a broader transition to flooding eigen- 
states. 



VII. SUMMARY AND DISCUSSION 

We provide a detailed quantitative description of the 
flooding of regular islands. By using the projection of 
eigenstates onto regular basis states, which defines the 
weights W m , the process of flooding can be described sep- 
arately for each torus. The distribution of these weights 
in the kicked system agrees accurately with the distri- 
bution obtained by the proposed random matrix model. 
This model depends on two parameters only: the cou- 
pling strength v between regular and chaotic basis states 
and the ratio of the number of those states A^eg/A^h. The 
connection of this coupling strength with the parameters 
of the kicked system is given by Eq. (12U1) . 

From the random matrix model we gain the following 
general insights into the flooding of the m-th torus in 
terms of its tunneling rate 7 m and the Heisenberg time 
t„ rh : 



i) 7m < 



T H,ch 



All regular eigenstates on the m-th torus 



exist. None of the eigenstates predominantly extending 
over the chaotic region has substantially flooded the m-th 
torus. 

ii) 7m = (27r) 2 1 : No regular eigenstates on the m-th 
torus exist. Some of the eigenstates predominantly ex- 
tending over the chaotic region have substantially flooded 
the m-th torus. 

iii) 7 m > Cj^j nT7 : ^ u °^ * ne eigenstates predom- 
inantly extending over the chaotic region have substan- 
tially flooded the m-th torus. 
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What do these results imply for the applicability of 
the semiclassical eigenfunction hypothesis? For a fixed 
system size in the semiclassical limit h Q g — > 0, which im- 
plies a roughly exponential decrease of j m , one ends up 
in regime i), in agreement with the semiclassical eigen- 
function hypothesis. In contrast, for small h e g 7^ fixed 
and systems with M cells and M — > 00, one obtains a 
large value for r H ch oc M, limited by dynamical localiza- 
tion only. Depending on the localization length one ends 
up in regime iii) for some or all tori m. As in our case 
one has r H ch ~ h e g/jo, regime iii) is realized for all tori, 
i.e. complete flooding of the island [ijj]. 

The universality in the transition from i) to ii) can be 
seen for the fraction of regular states f T eg(v) localized on 
a given torus. For the random matrix model this fraction 
does not depend on the ratio N reg /N c h and the agreement 
with the results for the kicked system is remarkably good 
for different quantized tori and values of h e g. In contrast 
to the disappearance of regular eigenstates on the m-th 
torus, the transition to flooding eigenstates on this torus 
is much broader and extends to regime iii). 

It is also important to discuss, what these results 
imply for the case of a single island in a chaotic sea 
[M = 1). Most commonly one is in regime i), i.e. regular 
and chaotic eigenstates exist and only mix at accidental 
avoided crossings. For a sufficiently small island, com- 
pared to the size of the chaotic region, regime ii) can be 
reached. Here h c g is small enough to quantum mechani- 
cally resolve the small regular island, but a corresponding 
regular state does not exist. It is not possible, however, to 
get into regime iii) where all eigenstates would be flood- 
ing eigenstates: In Eq. ([34]) we have h D g = 1/N and 
T H , ch = N c h w N such that the right hand side is approx- 
imately ir 2 N, which is always larger than the tunneling 
rates 7 m < 1. 

In the case of an island chain of period p embedded 
in a chaotic sea it might be possible to get into regime 



iii): In the derivation of Eq. (|34|) we now have to use 
v = N ch /(2N rcg ) m N/(2p) = l/(2ph e g), leading with 
Eqs. (H§ and iV ch « N = l/h e g to 7 > n 2 /(p 2 h c g). The 
right hand side can be smaller than 1 if p is sufficiently 
large while h c g is small enough to resolve the individual 
islands of the chain. Whether this is indeed possible in 
typical systems requires further investigations. 

This discussion shows that the semiclassical limit in 
generic systems with a mixed phase space, where islands 
of arbitrarily small size exist, is rather complicated. For 
example one can ask how small does h e g have to be such 
that at least one regular state exists on a small island 
of size A lcg ? Let us define the ratio r = h e g/A Teg The 
quantization condition Eq. (TTJ) implies that r < 2 is nec- 
essary to quantum mechanically resolve the island. How- 
ever, we find that the necessary ratio r becomes arbi- 
trarily small for small islands: Regime i) for m = 
requires 70 < l/r H ctl s» h c g. The tunneling rate 70 
is an approximately exponentially decreasin g fu nction 
70 ~ exp(— C/r) with C of the order of 1 flsL |32|| . Thus 
we have to fulfill exp(— C/r) /r < A reg , which for decreas- 
ing ^4 reg is only possible if r is sufficiently small. 

We conclude by emphasizing that the universality 
given by the random matrix model not only holds for the 
kicked system studied here, but is applicable to any sys- 
tem with a mixed phase space. The consequences for the 
semiclassical limit in the hierarchical phase-space struc- 
ture of generic systems needs much further investigation. 
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